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ABSTRACT 

A variational method that allows for replica-symmetry breaking is applied to directed 
polymers in an (A^ + l)-dimensional disordered medium. The noise studied here has 
gaussian correlations, i.e. it is short-ranged. In dimensions N < 2, the variational scheme 
yields only a strong-coupling phase and anomalous diffusion; while in dimensions N > 2 
it shows weak- and strong-coupling phases but no anomalous diffusion. Comparisons are 
made with the results of Mezard and Parisi [11] for noise with power-law correlations. 
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The sluggish rate of progress on strongly disordered systems attests to the difficulty 
of the problem. Some headway has been made over the past several years in the area 
of manifolds subject to a random potential. These may represent the most tractable 
models in which disorder plays a crucial role. Furthermore, they have been tied to various 
physical phenomena. For instance, the one-dimensional incarnation, a directed polymer 
in a random environment [1-12], has been related to domain walls in two-dimensional 
disordered ferromagnets [1], surface growth [13], and randomly stirred fluids [14]. 

By a "directed" polymer, one means a polymer or random walker that always proceeds 
in a positive direction along a given coordinate in an (A'"+ l)-dimensional space. The effect 
of the disorder is manifested by the walker's accumulating a series of random weights, one 
associated with each site through which it passes. The walker's attempt to maximize the 
benefit of the random potential results in far greater stretching along the axes transverse 
to the directed axis than the purely entropic spreading of its non-random counterpart. In 
the non-random case, the transverse fluctuations are diffusive, i.e. scale as: {u'^{£)) ~ £ 
(where £ measures the distance along the directed axis and uj that along the transverse 
axes); whereas in the random case, the transverse fluctuations are super diffusive, i.e. scale 
as {u'^{i)) ~ (for large i) where > ^ (where f{x) indicates averaging f{x) over 
realizations of the disorder). 

This nontrivial behavior is deeply rooted in the quenched nature of the disorder — 
as is the difficulty in the analysis thereof. Thus, it seems natural that the issues and 
techniques familiar from another problem featuring quenched randomness, namely spin 
glasses, have resurfaced here. Derrida and Spohn [6] made the connection to spin glasses 
explicit when they uncovered a mapping between directed polymers on a Cayley tree and 
the Random Energy Model, which is in some sense "the simplest spin glass." [14] 

The replica approach, another technique from the study of spin glasses, has made its 
way into the study of manifolds in disordered media. Early on, Kardar [4] used it to suggest 
that {lij^{£)) ~ £^ for (1 -|- 1) -dimensional directed polymers with delta-function correlated 
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disorder. More recently, Mezard and Parisi (MP) [11] applied a variational version of the 
replica method to the general problem of random manifolds pinned by quenched random 
impurities. This variational approach has also found applications in the study of (non- 
directed) polymers in a random media [16], protein-folding [17], disordered vortex lattices 
[18], the random-field Ising model [19], and interference effects in variable-range hopping 
[20,21]. 

In their extensive treatment of manifolds [11], MP included the case of directed poly- 
mers with delta-function correlations by invoking dimensional analysis, that is, by consid- 
ering power-law correlations that scale in the same way as the delta functions: for instance, 
l/|a;| would replace 5^^\x). In the present work, we apply the variational replica method 
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to directed polymers with noise correlations that fall off like a gaussian ~ e~'^ — a genuine 
short-range behavior. 

To recap briefly, MP [11] studied a very generic hamiltonian for D-dimensional man- 
ifolds in an {N + D)-dimensional space: 

H[uj] = IJdxY, (^) +Jdx Vix,ujix)), (1) 

where a; (a;) is an A'"-component vector field and a; is a D-dimensional vector. The potential 
V{x,u) is a random variable with zero mean and correlation: 

Vix,u)Vix',u') = -iV5(^)(a; - x')f (^^^^9^) " 

The scaling with N was introduced to facilitate the large- A?" limit in which the variational 
approach used is expected to become exact; it will be maintained here for purposes of 
comparison. 

MP considered noise correlations of the power-law form: 

(^) ~ 2(1^ (^) ' 
for uj^ ^ 1. Under these conditions, their calculations split quite naturally into two cases: 
noise with "long-range" correlation (either D > 2 or D < 2 and 7 < ^^"^ noise with 
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short-range correlation {D < 2 and 7 > ^z^)- Again, we note that MP only considered 
directed polymers with delta-function correlations indirectly by relating them to power-law 
correlations with the same scaling, in particular, by setting 7 = l + N/2. One might find it 
somewhat dissatisfying that this prescription places directed polymers with delta-function 
correlations in (1-1-1) and (2-1-1) dimensions into their "long-range" category. Another 
possible source of discontentment is MP's treatment of "short-range" case, in which they 
expanded the short-range behavior of the correlation in a power series, truncating it after 
only two terms. The present calculation employs a bona fide short-range form for the 
correlation and requires no such truncation. It will be shown to reproduce the main 
results of ref. [11] for directed polymers. 

In the directed-polymer case {D = 1), the partition function corresponding to H[u] 
(eq. (1)) is a path integral, 

»{aj,x} 
'{0,0} 

(where /3 is the reciprocal temperature) which can be shown to obey the equation: 

dZ{u!, x) 



/ Vu expi -(3H[u;] \ (4a) 



dx 

where Ti. is the (random) operator: 



= -n{u;) Z{u;,x), (46) 
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= - + PV{x,u{^)). (4c) 

i=l * 

But it can be usual first to replicate the system n times and average over the disorder; one 
then obtains a path integral which obeys 



az-(K},x)) 



= -^n(K}) Z^{{uja},x), (5a) 



dx 

with the following (non-random) n-body hamiltonian 

i=l a=l a^b ^ 



t^a - t^b)^ ) , {5b) 



with attractive interactions given by / (as defined in eq. (2)). 

The interesting exact solutions to such an n-body Schrodinger-hke equation have 
aheady been exploited, [4,8] and standard perturbative approaches have been shown to 
yield to unphysical results. [11] Hence, we shall proceed here with a variational scheme 
of the Rayleigh-Ritz type. That is, we shall study the ground-state energy of Hn{{(^a}) 
using the following variational wave function: 



This method of searching for the variational ground state of an operator hamiltonian is 
convenient when considering one-dimensional manifolds {D = 1) [12,17,20,21] and differs 
slightly from the approach taken in ref. [11] where higher dimensional manifolds are 
considered as well. The results, however, do not differ. 

Keeping in mind the eventual n — > limit and the possible emergence of replica- 
symmetry breaking (RSB) therein, a gaussian form was chosen for calculational ease prior 
to that limit. The matrix m of variational parameters is chosen as an n x n hierarchical 
matrix with i^-step breaking (as introduced by Parisi [22]). In the one-step {K = 1) 
breaking scenario, one can envision breaking the replicas into groups with the couplings 
between "particles" belonging to the same group differing from the couplings between 
"particles" belonging to different groups. Next, two-step breaking divides the groups into 
subgroups, and so on, leading to the hierarchical (ultrametric) structure. Eventually, the 
limits — > oo and n — > are to be taken. In addition to offering one more variational 
parameters (than the replica-symmetric version), it is hoped that the variational scheme 
with RSB can better mimic situations in which two or more diverging paths have nearly 
degenerate energies, including those with the very lowest energies. [9] 

With the gaussian wave function, calculating (Tin)* requires computing Trim) (which 
arises from the expectation of the kinetic-energy portion of Tin) and certain functions of 
(m~^ + "^6"6 ~ ^"^a I) (which emerge from the interaction terms). To be more specific, a 




n 



N 



(6) 



a,b i=l 
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Taylor expansion and Wick's theorem applied to an interaction term lead to [11]: 



N 



j\NJ 



fU) (Ar + 



2j - 2)!! (™a,a + ^6,6 " 2™a,b)" 



/ 



j\Nj {N - 2)!! 



2J 



(7a) 



where 
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do; q: 2 ^e / 



r(f ) Jo ' \ N 

We are currently interested in the case in which the correlations have a gaussian form 



(76) 



N 



[1 



N 

2 \ — 



ttN 



exp 



N 



(8a) 



where N^/ij, serves as a correlation length and fg{z) is normalized here so that it ap- 
proaches a delta function in the fj, —>■ oo limit. The corresponding / is: 



N 

2 \ — 



TVN 



1 + 



N 



N_ 
2 



(86) 



For power-law correlations (eq. (3)), the same procedure yields: 

9r(i + f-7)/2.^'-^ 



fp.ti^) = 



2(i-7)r(f 



iz 

N 



(9) 



provided 7 < 1 + Note that fg{z) and fp.i.{^) have the same large-z behavior when 
7 = 1 + Yi furnishing evidence for the validity of MP's dimensional-analysis approach 
within this scheme. (Setting 7 = 1 -|- ^ requires some small- 2; regularization; otherwise, 
the coefficient of fpj, would diverge.) 

The K-step hierarchical matrix rh has K + 2 parameters {a, ao, ai, . . . , ax} which take 
on the form [a, a(tt)], where a{u) is a function on the interval [n, 1], in the n — > limit. 
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[22] Ref. [21] showed how the eigenvalues of m, which take on a similar form [A, A(tt)], are 
related to [d,a{u)]. In terms of the eigenvalues, {'Hn)<i! becomes: 



Nn 

where 



Q(u) = r^X-\v) + -X-\u). (106) 
Ji f u 

The first two terms in eq. (10a) correspond to Tr(m), i.e. the sum over the eigenvalues. 
Note that du/u^ is roughly the degeneracy of X{u). 

Now comes the variation. One obtains solutions for the "best" A(u) through the 
following procedure. First, take a functional derivative of eq. (10a) with respect to A(tt) 
and set it equal to zero; this yields the stationarity equation: 

X\u) = 4(3'yyvf'[Q{v)]+uf'[Q{u)]'^, (11) 

where /' = df/dQ. As a step toward finding an equation which is local in u {i.e. no 
integrals over u), take a derivative with respect to u, which results in: 

2X{u)Xu{u) = -4p^f"[Q{u)]X-\u)Xu{u), (12a) 

where A^ = dX/du. This result implies that either Xu{u) = or 

2-^r')^Hu) = -f"[Q{u)]. (126) 

To pursue the latter case, put in the desired form of f"{x): 



1 + ^ 



Af+4 



(13a) 



and for simplicity let 

N KttN 



(136) 
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and 



t = (13c) 



N+A 
2 



These substitutions produce: 

2-'^p-^\^{u) = s\l + tQ{u) ' (14a) 
which upon inversion becomes: 



2 

JV+4 



2-'p-'s-'X%u) = l + tQ{u), (146) 

where Q{u), recall, is given by eq. (10b). Now, eq. (14b) is still non-local in u, so take 
another derivative with respect to u, furnishing: 

" {nTi) (2/^'^)"^ A-^(«) A.(^.) = - U-\u) X^{u). (15) 
The solution of eq. (15) is either Xu{u) — or: 

JV+4 

Note that eq. (16) and X{u) = Const, (which corresponds to Xu{u) — 0) are merely the 
possible local solutions. One must piece together a function on the entire interval [n, 1], 
consisting of these local solutions, which satisfies the (non-local) stationarity equation (eq. 

(11))- 

The pattern of RSB found by MP in their "long-range" solution fuses together a 
power-law behavior (similar to eq. (16)) for n < u < Uc and a constant for ttc < tt < 1. 
MP demonstrated that this small-w, power-law dependence is intimately connected to the 



property of superdiffusion: (a;^^^^)^ ^ £2;/ -^^^j^ ^ > ^- In particular, they found for ly the 
"Flory" exponents [5] {e.g. = | for directed polymers in 1 + 1 dimensions, instead of the 
exact result v = | [3,4]). 

For the gaussian problem, we will search for solutions of a similar nature. Note, 
however, that for N > 2, the power-law is not well behaved for small u in the n — > limit 



— then we will seek solutions like those derived by MP for their short-range case; which 
consist of two constants {i.e. one-step breaking). Hence, the present calculation also splits 
into two parts: N < 2 where one finds the full- breaking, power-law solution for small u 
and which displays anomalous diflFusion; and N > 2 where one finds the one-step breaking 
solution, a phase transition (see below), but no anomalous diffusion. 
First, for the case N < 2, let us assume a solution of the form: 



Xiu) 



{2p^sy 
{2phy 



6u 
{N+4)t 



iV+4 
2-JV 



iV+4 
2-JV 



if n < u < Uc', 
if Uc < u < 1; 



(17) 



then one can determine Uc by considering eq. (14) at tt = 1 (bypassing some tedious 
integration). This procedure yields the following equation for Uc- 



N 

~2 



1^ 



N{N + 2)l3" 

'a ^ 



6Ur 



N + 4 



+ 



6Ur 



(iV + 4) 



- 1 



0. 



(18) 



This solution, eqs. (17) and (18), was also obtained numerically by solving iteratively a dis- 
cretized version of the stationarity equation (eq. (11)), giving credence to the assumptions 
made in (17). 

Note that here Uc is /^-dependent. This feature differs from MP's solution for power- 
law correlations, where Uc is a constant. [11] With a /9-dependent Uc, one ought to check 
for solutions with Uc lying outside the interval [0, 1] , as that might indicate some kind of 
a phase transition, but Uc always lies within this range for A/" < 2. When = or = oo 
(the delta-function limit), the first term in eq. (18) is zero, and one finds Uc = "^^g-^, 
which is precisely the constant found by MP for power-law correlations (if 7 = 1 -|- N/2). 
Then, as /? increases (or decreases) Uc approaches though does not reach zero, so that 
< u < For finite /j, and large one finds Uc ~ This property distinguishes 

predictions concerning the large P behavior of quantities (such as the free energy) made 
by the gaussian and power-law calculations. 
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For N > 2, we look for a solution with one-step breaking of the form: 

X / _ / 0, if < w < Mc; 
^^""^ ~ [A, if«c<«<0. 

Substituting this solution into the hamiltonian (eq. (10a)) yields: 



(19) 



Note that when X{u) ~ 0, Q{u) — oo. Next, taking derivatives with respect to A and Uc 
and setting them equal to zero results in: 

JV JV+2 

and 

JV JV 

the stationarity equations. 
Eliminating A leads to: 

Again, we check for situations in which Uc lies outside the interval [n, 1] . This time we 
find some. In fact, we find /?c, the critical simply by setting ttg = 1 in eq. (23): 



t3 /ttN 

Y V7 



JV 

2\ - 
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/^c = -(-^) {N-2) . (24) 



For high temperatures, (3 < I3c, the solution of eq. (23) has Uc > 1; so then the solution is 
X{u) = for the entire interval [n, 1] — the (trivial) replica-symmetric solution. Here too, 
the assumptions made in eq. (19) were tested using the discretized version of eq. (11). 

The presence of a transition to a weak-coupling phase for iV > 2 is in accord with the 
known behavior of directed polymers. [6] As mentioned earlier, in their consideration of 
short-range correlations, MP expanded f{z) fa fo + fiz + ■ ■ ■ for small z, keeping only the 
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first two terms in their expficit calculations. In the gaussian case, no such approximation 
was deemed necessary, and yet the same basic results were obtained. Hence, it indicates 
that MP's results are not purely an artifact of the truncation procedure. 

In summary, we have applied a variational replica method to (N + 1) -dimensional 
directed polymers with gaussian-correlated noise. The calculation divided quite naturally 
into two regimes. In the first, N < 2, the variational solution displayed full replica- 
symmetry breaking with a power-law behavior for small u. This feature is associated with 
superdiffusion, and in particular the "Flory" exponents for the transverse fluctuations 
are found. Similar behavior was observed by MP [11] in their so-called "long-range" 
solution for noise with power-law correlations. Another trait shared by the gaussian- 
noise and power-law- noise calculations [11] is the lack of a transition to a weak-coupling 
(high-temperature) phase for N < 2. On the other hand, one distinction between the 
two is their dependence on the inverse temperature especially for large In the 
second regime, N > 2, the variational solution displayed a phase with no replica-symmetry 
breaking for high temperatures and a phase with one-step replica-symmetry breaking at 
lower temperatures, but no anomalous difi^usion. 

The variational scheme manages to uncover the phase transition where it is expected 
(AT > 2) and the absence of a transition where it is expected {N < 2); that is, it obtains 
the correct phase diagram. It predicts anomalous diffusion for N < 2; however, it fails 
to flnd the correct exponent u for the transverse fluctuations. Moreover, it predicts no 
superdiffusion in the strong-coupling phase for A?" > 2 as is anticipated from the simulations. 
[10] One scheme for furthering calculations involves an expansion (in 1/A^) around the 
solution at N = oo which is proposed to be exact. [23,24] Another conceivable (albeit less 
controlled) approach would be to improve the gaussian assumption for the variational wave 
function, perhaps by using a Lanczos-type approach [24] to go beyond the basic Rayleigh- 
Ritz method. It is hoped that the present work has resolved a few of the nagging issues 
found in MP. However, the more important issues remain unresolved and more work is 
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required before this chapter can be cfosed. 
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